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Abstract 

In this work we study in a formal way the density dependent hadron field the- 
ory at finite temperature for nuclear matter. The thermodynamical potential and 
related quantities, as energy density and pressure are derived in two different ways. 
We first obtain the thermodynamical potential from the grand partition function, 
where the Hamiltonian depends on the density operator and is truncated at first 
order. We then reobtain the thermodynamical potential by calculating explicitly 
the energy density in a Thomas-Fermi approximation and considering the entropy 
of a fermi gas. The distribution functions for particles and antiparticles are the 
output of the minimization of the thermodynamical potential. It is shown that in 
the mean field theory the thermodynamical consistency is achieved. The connection 
with effective chiral lagrangians with Brown- Rho scaling is discussed. 



PACS number(s): 21.65. +f, 21.30.-x, 25.70.-z 



1 Introduction 

Recently effective field theories have been often used in order to describe hadron proper- 
ties since the solutions for the fundamental theory of strong interactions, the Quantum 
Chromo Dynamics (QCD), are still an unreachable problem at the hadron physics energy 
scale. The study of nuclear matter and finite nuclei properties at finite temperature has 
become an important problem since a large variety of data, where matter is being tested at 
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extreme conditions of density, pressure and non-zero temperature, are becoming available 
in the modern experimental facilities which are already operational. 

The motivation for the formulation of Density Dependent Hadron Field Theory (DDHFT) 
m 12] may be linked with the relativistic Dirac-Brueckner Hartree-Fock (DBHF) theory 
Recently, one approach that has being used jUIS] is to consider the DBHF calculations 
for nuclear matter only as a guide for a suitable parametrization of density dependence 
of the meson-baryon coupling operators. In fact, the parameters are adjusted in order 
to describe the nuclear matter and some finite nuclei properties, as it is done in the 
usual Walecka model parametrizations [U]. In the present work we discuss how to treat 
the nuclear matter in the DDHFT at finite temperature and show that, in the mean 
field approximation (MFA), the thermodynamical consistency is achieved. This is in line 
with previous calculations done on the same subject, where thermodynamical consistency 
was assumed fj^. We also show that the mean field approximation is equivalent to the 
Thomas- Fermi approximation (TFA). We follow two different methods in obtaining our 
results. We first start from the calculation of the grand partition function in terms of the 
density dependent Hamiltonian, which is truncated at first order. Once the thermody- 
namical potential is found, all other thermodynamical quantities can be obtained from it. 
In the second approach, we start by writing the thermodynamical potential in terms of 
the energy density calculated in the Thomas-Fermi approximation and the entropy of a 
fermi gas. The distribution functions for particles and antiparticles are the result of the 
minimization of the thermodynamical potential. 

Furthermore, there is a class of models, including those based on effective chiral la- 
grangians endowed with the Brown-Rho scaling, to which our formalism and our conclu- 
sions can be easily extended. This point is discussed later in this paper. 

The paper is organized as follows: in section 2 the DDHFT is presented. In section 3 a 
mean field approximation is performed and some termodynamic quantities are calculated. 
In section 4 the Thomas-Fermi approximation is used and the results are compared with 
the ones obtained in section 3. Finally, in the last section the conclusions are drawn. 
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2 Density Dependent Hadron Field Theory Formal- 
ism 



We start from the density dependent hadron field theory Lagrangian density of the well 
known Walecka-type models [HlllEl: 



7. - r.(p)^^ - ^T-r- I (1 + rs) A^) - (M - r^(p)0) 



C = ip 

-\p,u ■ r + \mlp, ■ r - \f,uF^'' , (1) 

where 0, tu^, and are the scalar-isoscalar, vector-isoscalar and vector-isovector 
meson fields and the photon field respectively, u^^ = d^uj^ — d^u^, F^y = d^A^ — d^A^, 
and = d^py — d^p^ — V p (p^ x p^)^ M is the nucleon mass, m^, m^, nip are the masses 
of the mesons and rs = — 1 (ra = 1) for neutrons (protons). The nucleon- meson density 
dependent coupling operators T^{p), r^(/)), Tp{p) are taken to be dependent on a Lorentz 
scalar functional, piip^if). In what follows we assume the vector density dependence 
description [T], i. e., (P = j'^j^, where jp = ip'jpip- Actually this has been the functional 
dependence most often used in recent applications [110 El- Notice that other possibilities 
for the functional dependencies of the coupling operators may be considered as well, for 
example, has been taken to be a function of ps = ipip in JQ. 

In this work our emphasis is to consider hadronic models in the context of the Walecka 
model. However, our formalism can be applied to a larger class of systems, i.e., any La- 
grangian density of the type £ = Cf + Cb + Cint, with Cp = ijj {i'y'^dp — M) ip correspond- 
ing to the fermion sector, Cb to the free bosons and Cmt being any interaction which can 
written as Cint = ^a{p)'4'^a4'V°', the summation running over the bosons, 6q, contain- 
ing Dirac and/or isospin matrices and V"' being any boson field (scalar, vector-isovector, 
vector-isoscalar, etc). 

As usual, the field equations of motion follow from the Euler-Lagrange equations. 
Here some care has to be taken since the coupling operators are now dependent on the 
baryon fields ip and ip through p. Hence, if C, stands for a generic field, the Euler-Lagrange 
equations are: 

dC ^ ( dC\ dCdp 



When the partial derivatives of L are performed relatively to the fields i\) and ■?/', they 
yield extra terms due to the functional dependence of the coupling operators in /)(■?/','?/'). 
These new terms are absent in the usual Quantum Hadrod Dynamic (QHD) model jH]. 
So, the equations of motion for the fields read: 

(9^9^ + m5)0 = r^(p)^^, (3) 
d^^^^'^mluj^ = TUp)^l^i^, (4) 



d^p'-' + mlp'' = ^J^i^T^^, (5) 



= |^(l + r3)7'^^, (6) 
[^,itd^-W{p))-M*]^P = 0, (7) 

where M* = M — T^{f))(l). Notice that in the equation of motion for the baryon field 
the vector self-energy consists of two terms, = S^^*) + S^, where: 

= r.(/3K + ^^r-p- +|(l + r3)A,, (8) 

where is the usual vector self-energy, pu^ = with = \ and, as a result of the 
derivative of the Lagrangian with respect to p a new term appears, S^, which is called 
rearrangement self-energy and will play an essential role in the applications of the the- 
ory. This term guarantees the thermodynamical consistency and the energy-momentum 
conservation, i.e., d^T^^^ = of the density dependent effective models, where the energy- 
momentum tensor is given by P]: 
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^^i-l^d"^ + + d^'^x^^^ + ■ P^^ + d^'AxF^^ . (10) 

The Hamiltonian operator follows from the T^o = ^ component of the energy-momentum 
tensor in the form: 



H = \ d^x Tn 



'00 



1 d^x (_a . (,v + S(p)) + /5(M - r<^(p)0) + So) ^ - V7^Sj^ + ^Fa^F^" 
(a'^c/.a.^ - mj^^) + 1.;,,.;^'^ - ^m>,^^ + \px, ■ ^ - \mlpx ■ . (11) 
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Formally, the thermodynamics of the DDHFT follows from the calculation of the grand 
partition function: 



Zq = Tr exp 



-/3 H{p) - ^ fiA 



i=p,n 



= exp [-pn] 



'12) 



where the trace goes over a complete set of states in the Fock space, (3 is the inverse of 
the temperature and Bp and Bn correspond to the proton and neutron number operator 
respectively, such that the baryon number operator is B = Bp + Bn- Although the 
equations of motion, cqs-lj^HZj), seem to be similar to the usual QHD ones [H], it should be 
noticed that the source terms for the meson fields contain in-medium correlations through 
the density dependence of the vertices and the rearrangement term contributes explicitly 
for the baryon field equation. Of course the solution of these field equations, as they 
stand, is a formidable task. For this reason, we are forced to consider approximations in 
order to obtain a tractable problem. Hence, we consider the DDHFT in the mean field 
approximation for nuclear matter at finite temperature. 



3 DDHFT for nuclear matter in the mean field ap- 
proximation 

We assume in this work, for the description of the nuclear matter, static meson fields, 
rotational and translational symmetry, charge conservation and no Coulomb field jU]. 
Therefore, all time and spatial derivatives vanish and only time-like components of the 
(J), uj and p meson fields remain. In the mean field approximation, the meson fields 
are replaced by classic condensed fields and the time-like components of the 0, oj and p 
condensed meson fields are called 0O; "-^o and respectively. In this approximation, only ip 
and if) are taken to be quantized field operators. The meson fields are obtained by thermal 
averages of the corresponding sources. The involved approximations are discussed later. 
In the nuclear rest frame, consistently with the MFA, we have p = ip'^o'^- Under these 
approximations the Hamiltonian of eq. (fTT|) reads: 

H = Jd^xH 

= J d'x {^t . V + /?(M - r^(p)0o) + ^ 

+lml<Pl - ^mlul - ^mlbl^ (13) 
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where 



sr = r^(p)^o + 



rp(p) 



(14) 
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It is important to emphasize that in the MFA the rearrangement term cancels out in the 
energy density, as a consequence of the cancellation of this term in the exact Hamiltonian 
of equation (jllj) . since the density dependent interaction terms in the Lagrangian density 
of equation (Q) involves no derivatives with respect to p. 

Now we proceed to calculate the partition function. In order to perform the trace 
calculation in (fT^. we follow the same procedure used in the calculation of quantum gases 
with density dependent interactions and in the molecular field approximation to the 
Heisenberg Hamiltonian JT]. For the evaluation of the trace in eq. fll2|) it is necessary 
that, in a convenient basis of states, the exponential function decomposes itself in a sum 
of independent terms. As it is seen next, this can be done only when we approximate 
the Hamiltonian operator in eq. (jl3p by a one-body operator, obtaining a linear term in 
the baryon number operator. Analogously to ordinary statistical mechanics problems jTUI 
we take the Hamiltonian expanded around the equilibrium mean density po = (p)- We 
then keep only terms up to the first order in p — Pq. Higher order terms are going to be 
neglected. 

In the DDHFT it is assumed that the functional form for the density dependent 
coupling operators is sufficiently well behaved as a function of p. In eqs. P3II14|) we expand 
the vertices around po in the form. 



and substitute the normal ordered scalar and isovector density by their thermal averages 




or,. 



(15) 



ijjljj 



(ipij) = ps , 

(^^370^) = P3 



In this way we obtain, for the Hamiltonian density up to first order in (p — po)- 




(16) 



where 




1 dT,{p) 

2 dp 



bop3 - 



dp 



P=PO 



0OP5 
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To calculate the trace in eq. (jl2|) . we use the single particle basis associated with the Dirac 
field equation, eq.Q, calculated at density po, where po = B/V is the equilibrium nuclear 
matter density. This equation can be solved exactly if one seeks stationary solutions of 
the plane wave form, ip = exp(zp-x — z£:(p)t), where ip{p) is a four component spinor. 
Such procedure is well known in the literature and details can be found in references jUJ E] . 
So, in the MFA, the baryon spinors are eigenvectors of the stationary Dirac equation: 

i-ta • V + /?(M - r^(po)0o) + So(po)) ^ = eilj, (17) 

and the baryon energy is given by: 

e{p) ^ 5^ = T^uo + ^T,bo + ± (p2 + M*y^' , 

where M* = M— r<^(po)0o- Going through the usual steps to solve the Dirac equation, one 
finds that e^^\e^~^) corresponds to the baryon m(p) (anti-baryon f (p)) four component 
spinor. Proceeding as in the ordinary Walecka model ^J, both positive and negative 
energy states u and v are found and the field operator, in the Schrodinger picture, can be 
expanded as: 

V^ = ^(n(p,A)aJ,, + t;(p,A)6y . 

pA 

In the summation A stands for the spin and isospin projections. The operators a^^^ 
and bp)^ are interpreted as creation operators for baryons and anti-baryons and satisfy 
anticommutation relations 0. Thus, substituting the field operators into eq. (fT^ and 
taking into account the expansion given in eq. (|16|) the normal ordered Hamiltonian takes 
the form: 

H = ^(p2 + M*2)'/Va;,ap, + 6j,,6p,) + (r^a;o + S«(po))s + ^6o4 

pA ^ 

+V (^ml<Pl - ]^mlul - ^my, - So^(po)po) , (18) 

where B:^ = Bp — Bn, V is the volume of the system and the nucleon number operator is 
given by 

ps 

Since in the present approximation all operators in the exponential function defining the 
partition function, eq. p2p . are diagonal in the basis of eigenstates of the baryon and 
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antibaryon number operators, the grand partition function can be exactly calculated. 
The results are analogous to those for the Walecka model [U] reading 

Q{(3,fip, fin, (po,^o, bo) = —InZc 

V 2_ 
Y ■ 

I %=p,n P 

+ ln(l + e-^(^*(P)+'^*))] , (19) 

1 /2 

where E*{p) = (^p^ + M*^^ and the effective chemical potentials, n* and /x* , are defined 
as 

r 

* T-\ P u 

f^p — f^P ~ ^ u}^0 — Zjq , 

z^: = /i„ - r^o^o + - (20) 

Notice the explicit contribution of the rearrangement term in the above expressions. The 
pressure is obtained through the relation P = —Q/V and the energy density by: 

^ E ld{(3n) 

t = — = — + flppp + /inPn ■ (21) 

The finite temperature meson field equations can be determined by assuming that the 
thermodynamic potential VL be stationary for variations of these fields. This is to be 
expected for a system in thermodynamical equilibrium. So, 

dVL 



yield the coupled equations 



0, ^ = u;o,&o,0o 



m|0o - r<^(po)Ps = 0, (22) 
mlujQ - r<^(po)Po = 0, (23) 
-^>o-^^P3 = 0, (24) 



where the thermal scalar and baryonic densities are defined as 

d^p M* 



= 2 (/.. + /.-). (25) 



Ps 

i=p,n ' 

Po = (^tV) = Pp + Pn , (26) 



(27r)3E*(p) 

'n , 

P3 = (V'7°^3^) = Pp- Pn , (27) 
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with the distribution functions given by 



fi± = r/ ^ / N , , i = p,n (2J 



and the proton and neutron densities pp and p„ reading 



A = 2 / ^(/.+ -/.-), ^=p,n. (29) 



The equations of motion, eqs. ()22ll!Mjl . are completely consistent with the MFA. They 
follow from eqs.(jniIZ|) when the sources are approximated by their thermal averages as 
discussed above. For example, for the scalar field source, T^{p)ijjip r(^(po)p5, etc. 

Next, we reobtain, within the Thomas- Fermi approximation, the same expressions just 
obtained for the pressure, energy density and meson fields, demonstrating the consistency 
of our approximations. 



4 DDHFT in the Thomas- Fermi approximation 

We first define the thermodynamical potential, following the notation in ^2), as 

^ = S - TS - PpPp - PnPn , (30) 

where S are the energy and entropy density respectively, T is the temperature, Pp (pn) 
is the proton (neutron) chemical potential and Pp and p„ are respectively the proton and 
neutron densities as given in eq. lj^ . calculated in such a way that p = pp + pn- Here, the 
distribution functions /j+ and /j_ for particles and anti-particles have to be derived in 
order to make the thermodynamic potential stationary for a system in equilibrium. The 
entropy density S is obtained from 

Next, we obtain the energy density in the TFA, where the conserved energy-momentum 
tensor, eq. lfTmi . becomes: 



(32) 



^From the energy-momentum tensor one easily obtains the Hamiltonian operator, which 
is the same as already displayed in equation (fTSj) and from it, the energy density operator. 
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which reads [U]: 

pA ^ pA 



psra 



where A = {s, t^} are the spin and isospin indexes. In obtaining the above equation we 
have assumed that the expectation values of the meson fields are constant classical fields. 
Moreover the Hamiltonian density is simply T-C{po), contrary to the expression (|T6|l used in 
the calculation performed in the last section. The energy density in the TEA is obtained 
through the substitution of the baryon and antibaryon number operators by their thermal 
averages. So, the energy density can be written in the semi-classical TFA as 



i=p,n 



+ ^00 - - -Y^o ■ (34) 

After straightforward substitutions, eq.flHUI) becomes 



i=p,n 



-■^T. (35) 

i=p,n \ ' 

For a complete demonstration of the above shown expressions obtained in a Thomas- 
Fermi approximation for the non-linear Walecka model, please refer to At this point, 
eq. ()35|) is minimized in terms of the distribution functions for fixed meson fields, i.e., 

=0 z^j. (36) 

fi~Jj±,<t>0,^0,bo 



5f. 

For the proton distribution function, the above calculation yields 



E*{p) + - + r^a;o + ^60 = -T In (^^^j • (37) 

Similar equations, with some sign differences are obtained for the anti-proton, neutron 
and anti-neutron distribution functions. The effective chemical potentials are the same as 
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defined in eq. (|2Up and from tlie minimization we reobtain for the distribution functions 
the expression given in eq. (j28|) . In other words, with the distribution functions given 
in substituted into equation (jH^j) . equation lfT^ is exactly reproduced. In the above 
calculation we have used po and ps as defined in eqs.(|^B|) and respectively. Within 
the Thomas-Fermi approach the pressure becomes 



i=p,n y p 

Again we recover our previous expressions for the pressure and the energy density. It is also 
important to stress that the thermodynamical consistency, discussed in pQ, which requires 
the equality of the pressure calculated from the thermodynamical definition i7 = —PV^ 
where f2 is given in equation (jl9p or (j35|) and from the energy-momentum tensor, is also 
obeyed by the temperature dependent DDHFT. This can be easily verified calculating the 
pressure from the thermal average of the energy-momentum tensor ()32|) in the MFA as 



P = \jlk%^- (39) 



5 Conclusions 



In this paper we have incorporated temperature effects in the DDHFT formalism in- 
troduced in [21 E] . Two different approximations were performed, namely, a mean field 
approximation where the coupling operators are expanded in powers of the density and 
the Thomas- Fermi approximation. In the first case, the Hamiltonian is expanded in pow- 
ers of (p — Po) and the Dirac equation is solved yielding the baryon energies. The grand 
partition function is calculated and the thermodynamic quantities are derived. In the sec- 
ond case, the thermodynamical potential is the starting point, where the energy density 
is calculated in the semi-classical TFA. The thermodynamical potential is minimized in 
terms of the unknown distribution functions, which are then obtained. Thermodynamic 
consistency is proved. 

As a final remark we discuss briefly another possible application of the formalism we 
have derived above, i. e., the study of effective chiral lagrangians with Brown- Rho (BR) 
scaling. In 1991 Brown and Rho (BR) proposed an in-medium scaling law for the 
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masses and coupling constants for effective chiral lagrangians. In ref. ^3] the authors 
proposed an effective Lagrangian whose parameters scale in nuclear medium according to 
the BR scaling. In its simplest version the Lagrangian in the MFA is |14t [T3] 

£ = i^b,{id^-g:{p)uj^)-M*{p) + hcP]^lj 

where in the notation of ref-jTH] is the nucleon field, uj^ the isoscalar vector field, 
an isoscalar scalar field and the masses with asterisk are BR-scaled as introduced in 
ref.pni- The scaling of the vector coupling constant is left arbitrary and h is taken 
constant. So, the Lagrangian in ()40|) is of the form of a Walecka-like Lagrangian and all the 
finite temperature formalism that we have developed for the DDHFT can be immediately 
applied to these lagrangians. The thermodynamics of effective lagrangians with BR scaling 
has been studied in for zero temperature. The study of the validity of the BR scaling 
hypothesis for the non-zero temperature case is currently under investigation. 
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